Entangled Dynamics of a Stiff Polymer 
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Entangled networks of stiff biopolymers exhibit complex dynamic response, emerging from the topological 
constraints that neighboring filaments impose upon each other. We propose a class of reference models for 
entanglement dynamics of stiff polymers and provide a quantitative foundation of the tube concept for stiff 
polymers. For an infinitely thin needle exploring a planar course of point obstacles, we have performed large- 
scale computer simulations proving the conjectured scaling relations from the fast transverse equilibration to 
the slowest process of orientational relaxation. We determine the rotational diffusion coefficient of the tracer, 
its angular confinement, the tube diameter and the orientational correlation functions. 
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Semidilute solutions of stiff biopolymers form entan- 
gled networks with remarkable mechanical properties 
and complex dynamic response |Q]] ; examples include F- 
actin SEBJaSill, microtubules H, the fd virus ||, 
and xanthan IllOn . Such a behavior originates from topo- 
logical constraints imposed by the impenetrable neigh- 
boring filaments. A key to the dynamics of the individ- 
ual filaments is the reptation concept pioneered by Ed- 
ward s fllll and de Gennes fl2rj . and later extended to 
rods 11311 : it summarizes the complex interaction of a 
single polymer with its surroundings to an effective con- 
fining tube. Then, transport is restricted to sliding back 
and forth in the tube, which entails a continual remod- 
eling of the tube ends; the resulting snake-like motion 
was coined "reptation", from the Latin repere (to creep). 
Compared to the dynamics in dilute solutions, the relax- 
ation of the polymer from its initial position and config- 
uration becomes extremely slow. 

As a consequence, the dynamic processes of entan- 
gled solutions of, e.g., biopolymers cover many decades 
in time, posing a tremendous challenge both to experi- 
ments and simulations. For flexible polymers, the rep- 
tation concept is well established lfl4ll and fairly predic- 
tive fl5ll ; in the case of biopolymers, only the confining 
tube has been observed experimentally flrjl . Computer 
simulations of entangled polymer solutions encounter 
major difficulties to follow the reptation motion; yet 
they give insight into the relaxation within the tube 11711 . 
Progress beyond simple scaling arguments depends cru- 
cially on the design of generic models, which are as sim- 
ple as possible to follow the dynamics for sufficiently 
long times, yet complex enough to display key aspects 
of the underlying microscopic processes. 

In this Letter, we propose the following class of mod- 
els to explore single-filament transport in polymer net- 
works: consider the motion of a tagged polymer in a 
plane, entangled in a course of immobilized obstacles; 
the latter represent the topological constraints due to the 
neighboring filaments. The reduction of dimensionality 
still captures the physics of entanglement since the repta- 



tion motion is essentially one-dimensional lflll. [l2ll. Stiff 
polymers are rather straight and therefore can be em- 
bedded in a plane, neglecting the torsion of their space 
curve. The orientation of the confining tube is persistent 
on the longest time scale of interest; thus in video mi- 
croscopy experiments, the non-trivial reptation motion 
of a labelled polymer found initially in the focal plane 
takes place in this plane. As a benefit of the simplifica- 
tion, the computational complexity is lowered substan- 
tially, permitting a thorough investigation of slow dy- 
namic processes. 

This class of entanglement models sets a framework 
for the plethora of polymer aspects: depending on the 
specific problem at hand, various polymer models may 
be employed ranging from a flexible chain of beads and 
springs to a rigid rod. It is essential to characterize and 
understand several limiting cases. The physics of entan- 
glement is singled out in the limit of hard-core interac- 
tion, vanishing width of the polymers, and zero exten- 
sion of the obstacles. Then, all ramifications of excluded 
volume are eliminated, all configurations are permitted 
and equally likely, and all non-trivial dynamic correla- 
tions are due to entanglement. In particular, the limit 
circumvents the nematic phase transition. 



In de Gennes 's seminal paper on reptation 11211 . a spe- 
cific realization of this fundamental limit has been con- 
sidered: an infinitely thin, flexible polymer moving be- 
tween fixed point obstacles. Research in the last decades 
rendered this model a hallmark in reptation theory. Com- 
plementary to a flexible polymer with respect to the 
bending stiffness is a needle, i.e., a straight, rigid object 
of negligible width, characterized solely by its length L. 
In the remaining part, we focus on this important refer- 
ence system. In particular, we validate in detail the tube 
concept for rods lfl3ll and extend it towards a complete 
theoretical description of the rotational dynamics. 

For a three-dimensional suspension of needles of num- 
ber density n, Doi and Edwards lfl3ll conjectured an 
asymptotic suppression of the rotational diffusion co- 
efficient, D rot ~ n~ 2 as n — > oo. The slow dynam- 
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ics at long times is expected to be universal irrespective 
of the microscopic motion; so far, research focused on 
ballistic needles, i.e., without solvent. Early molecu- 
lar dynamics simulations of such needle liquids |[l8l. [l9ll 
show substantial deviations from Enskog theory. It was 
pointed out that the Doi-Edwards (DE) scaling of D mt 
is approached only very slowly and difficult to observe. 
Semiquantitative agreement was found using a pseudo- 
dynamics lEoll : yet, the DE scaling has not been val- 
idated by a simulation of the dynamics. Within an 
elaborate Boltzmann-Enskog theory l2lll . the onset of 
anisotropic diffusion from the dilute regime has been ex- 
plained recently. For fixed positions of the needles, the 
orientational degrees of freedom exhibit glassy dynam- 
ics mm. 

We have simulated the motion of a needle in a two- 
dimensional array of frozen, point-like, and hard obsta- 
cles; see Fig. [TJ The latter are distributed randomly, in- 
dependently, and uniformly in the plane with an average 
number density n. Then, the topology of the network of 
obstacles is characterized by the mesh size £ := n -1 / 2 , 
i.e., the typical distance between obstacles. The degrees 
of freedom of the needle encompass the center-of-mass 
position R and the unit vector of orientation u, the lat- 
ter being parametrized by a single angle 99. We com- 
pare ballistic and overdamped micro-dynamics of the 
needle. In the ballistic case, the total kinetic energy 
is conserved, and its value sets the overall time scale 
ro := L/v of the problem; v denotes the root mean- 
square velocity. For overdamped dynamics, the time 
scale ro is defined via the coefficient of unhindered dif- 
fusion along the axis of the needle, ro = L 2 /Z)jj; the 
free orientational diffusion is chosen in accord with hy- 
drodynamics, D° ot = 6D®/L 2 . In both cases, the model 
is specified by one single dimensionless control param- 
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FIG. 1 : Illustration of an entangled needle in a plane. The relevant 
length scales are the length of the needle L and the mesh size of 
the network £. The surrounding point obstacles confine the needle 
to a tube (shaded areas) of width d, a renewed tube (green) is tilted 
against the old one (red) by an angle e = d/L. 
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FIG. 2: Simulation results for the rotational diffusion coefficient, 
(a) Dashed lines are asymptotic fits to the predicted Doi-Edwards 
scaling, D^ t — A(n*)~ 2 ; the red solid line shows the result from 
a Boltzmann theory, D mt = 10.5/n*ro (24ll . (b) Deviation from 
the asymptotic behavior; note the different prefactors A in D^ t for 
ballistic and overdamped dynamics. 



eter, n* := nL 2 . Equivalently, the entanglement index 
pE := — log 10 (£/L) quantifies the relative importance 
of entanglement. 

Our simulations of the molecular dynamics are based 
on an event-driven al gori thm using a novel approach 
to collision detection II24I . I25I1 . Diffusion coefficients 
of the rotational motion have been extracted from the 
long-time behavior of the mean-square angular displace- 
ment (MSAD) 5ip 2 (t) := (A(p(t) 2 ) ~ 2D rot t and are 
shown on a double-logarithmic plot in Fig. |2^. In the 
investigated density range, the diffusion coefficient D mt 
varies over seven non-trivial decades. For dilute systems, 
n* <g 1, it depends on the micro-dynamics: In the ballis- 
tic case, the diffusion coefficient is suppressed in quan- 
titative agreement with a Boltzmann theory |[24ll . D mt = 
10. 5/ra* To. For overdamped motion, the needle is unaf- 
fected by the obstacles, and the diffusion coefficient is 
just given by D® oV Once the mesh size becomes compa- 
rable to the length of the needle, £ w L, the isotropy of 
rotational dynamics breaks down and a different trans- 
port mechanism develops |[2r3l . With growing entangle- 
ment, the needle is increasingly caged by the obstacle ar- 
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FIG. 3: Typical trajectories of the center of the needle for over- 
damped dynamics. With increasing entanglement, the needle is con- 
fined to a narrowing tube, and reptation dynamics emerges; see also 
the Supplementary Movie |25ll . 



ray, and eventually, its rotational motion is strongly hin- 
dered. Fig. [3] illustrates the emergence of reptation-like 
dynamics, accompanied by a drastic suppression of the 
diffusion coefficient D rot . For nL 2 > 10 2 , the data fol- 
low an asymptotic power law, D rot ~ (n*)~ 2 , over more 
than four decades in the diffusion coefficient. As a most 
sensitive test, Fig. |2j) compares D mt to its asymptotic 
behavior for increasing entanglement. Our results show 
that the mechanism of reptation is universal for ballis- 
tic and overdamped motion of the needle, which will be 
substantiated further in the subsequent analysis. 

The asymptotic suppression of D Iot is rationalized by 
employing the concept of a confining tube flUl . In this 
picture, the surrounding obstacles reduce the accessi- 
ble volume of the needle to a tube of diameter d and 
length L. The diameter is estimated from the require- 
ment to encounter no obstacles in the tube, d m 1 / nL = 
£ 2 /L. The constrained motion is illustrated in Fig.[3]and 
in a Supplementary Movie 12511 : the transversal and ro- 
tational degrees of freedom are essentially frozen, per- 
mitting only displacements along the axis of the tube. 
After traveling half its length, the needle is confined to 
a new tube tilted against the previous one on average by 
an angle e m d/L. The time ra to disengage from a 
current tube is estimated from the free longitudinal mo- 
tion inside the tube, t& ps to, independent of the density. 
Eventually, the orientation performs a random walk with 
step size e and constant rate l/rj; hence, the diffusion 
coefficient scales as 
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The given arguments apply likewise to ballistic and over- 
damped dynamics of the needle. In both cases, our data 
in Fig. |2] provide ample evidence for such a behavior, 
unprecedented in the literature. 

The quality of our data allows us to verify and quantify 
the assumptions of the tube model in detail on the basis 
of the MSAD and correlation functions; we will restrict 



the discussion to ballistic dynamics. The behavior of 
S(p 2 (t) is exhibited in Fig. from the very dilute up to 
the highly entangled regime with pE > 2. At short time 
scales, the motion is ballistic, 8<p 2 (t) = ( l f 2 ) t 2 for t < 

T CC) ti, with the mean collision rate r~j| = 0.845 n* /tq. In 
the dilute regime, the MSAD directly crosses over to dif- 
fusion. Entanglement effects emerge already at n* = 12 
(pE = 0.54): when t « r coll , the MSAD hits an interme- 
diate plateau, reflecting the angular confinement within 
the tube. Beyond this time scale, the transverse degrees 
of freedom are equilibrated. The MSAD increases fur- 
ther only after the tube is renewed at the time scale t&, 
and not until then, diffusion is observed. We use the 
square root of the measured plateau value as definition of 
e; it follows the expected scaling law over two decades, 
fixing also the pref actor, e = 1.3/n(L/2) 2 ; see Sup- 
plementary Fig. 1 125]. Nice data collapse is achieved 
for t > t co ii by rescaling the MSAD, (n*) 2 5ip 2 (t); see 
inset of Fig. @K. The disengagement time may quan- 
titatively be defined via extrapolation of the diffusive 
asymptote to the plateau, e 2 = 2D^t^\ our data yield 
D°& = 3.2/n 2 (L/2) 4 r , implying r d = 0.3r inde- 
pendent of the degree of entanglement. Supplementary 
Fig. 1 112511 also shows the tube diameter, inferred di- 
rectly from an analogous plateau in the transverse mean- 
square displacement; we find d = 1.3£ 2 /L and thus 
e = A.Od/L. 

More generally, the rotational motion in a plane is 
characterized in terms of orientational conelation func- 
tions Cfj,(t) := (cos|/tA^(i)]) for integer /j,. The caging 
inside the tube, 5tp 2 (t) = e 2 , should be reflected in a 
plateau close to unity, C^(t) ~ 1 — /i 2 e 2 /2 for interme- 
diate times r co i] < t < Td, which follows by a Taylor 
expansion of the cosine. For larger times, t > Td, one 
expects simple rotational diffusion; solving the diffusion 
equation on a circle yields C^t) ~ exp(— /x 2 D rot i). 
Thus, the initial orientation relaxes in two steps: fast 
equilibration to the plateau for t < r co u, and slow re- 
laxation from the plateau. We have checked that the 
measured C^t) indeed exhibit a plateau with the pre- 
dicted value, followed by an exponential relaxation with 
growing decay times r M = l/^i 2 D mt ~ {n*) 2 . Fig.Hh 
displays C\{t), the relaxation of the orientational per- 
sistence; this is the slowest process in the system corre- 
sponding to an equilibration time r rot := r\ = 1/D mt . 

Our results unambiguously prove the conjectured scal- 
ing relations for D mt , e, and d. Thus, the model reflects 
the generic DE scenario, demonstrating that the essential 
physics due to entanglement is captured. The predicted 
scaling behavior is, however, only observed in highly en- 
tangled systems with large entanglement index, pE > 1. 
In this regime, the trajectories indeed exhibit pronounced 
reptation with the typical sliding motion, see Fig. [3] 

For weaker entanglement, < pE < 1, the rota- 
tional dynamics is still suppressed due to topological 
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FIG. 4: Time dependence of the rotational motion. (a) The mean- 
square angular displacement of the needle develops a plateau at high 
densities, which scales as e 2 ~ (n*)~ 2 (cf. inset). (b) The expo- 
nential relaxation of the orientation correlation Ci (t) is preceded by 
a plateau very close to unity (see inset). Symbols indicate simulation 
results and solid lines fits to / exp(— t/r lot ); the fitted plateau value 
is consistent with the measured angular confinement, / = 1 — e 2 /2. 



constraints, but the DE scaling is obscured by crossover 
phenomena. The deviations from the predicted behav- 
ior are highlighted in Fig. |2j) by extracting the apparent 
amplitude of the power law. At pE = 1, where the fil- 
ament length already exceeds the mesh size by a factor 
10, the amplitude is still 15% below its true asymptotic 
value. One concludes that in order to observe the scaling 
with an accuracy of 1%, even stronger entanglement is 
required, pE > 1.6. 

In real polymer solutions, the entanglement con- 
straints are dynamically released and generated on a time 
scale comparable to t&. Within an extended model ac- 
counting for this renewal process of the obstacles, we 
have checked that the observed DE scaling is robust 12711 . 

For highly entangled networks, the finite filament 
width in experimental situations may become relevant. 
Eventually, a phase transition to a nematic order occurs 
for long rods at n^bL 2 ~ 1 as has been estimated by 
Onsager, where denotes the three-dimensional num- 
ber density of rods of diameter b. The density of obsta- 
cles of our two-dimensional representation is then calcu- 
lated to n « n-s&L, hence the nematic regime is expected 
for n*b/L = nbL > 1. For reconstituted F-actin solu- 
tions with filaments of L sa 50 um [16] and b = 7nm, 



we estimate that nematic effects are relevant only for 
n* > 7000 or pE > 1.9, provided one can neglect the 
small bending flexibility. 

A finite stiffness for the polymer introduces another 
length scale, the persistence length, quantifying the dis- 
tance over which the polymer appears as a straight rod. 
Due to thermal noise, there are transverse undulations 
which effectively blow up the width of the polymer. It 
is an open question if the finite flexibility assists for the 
tube remodeling resulting in an enhanced rotational dif- 
fusion, or if the additional effective volume leads to fur- 
ther slowing down. 
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